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Abstract 

We study the time evolution of entanglement of two spins in anisotropically coupled quantum 
dot interacting with the unpolarized nuclear spins environment. We assume that the exchange 
coupling strength in the z-direction J z is different from the lateral one J[. We observe that the 
entanglement decays as a result of the coupling to the nuclear environment and reaches a saturation 
value, which depends on the value of the exchange interaction difference J = \\Ji — J z \\ between the 
two spins and the strength of the applied external magnetic field. We find that the entanglement 
exhibits a critical behavior controlled by the competition between the exchange interaction J and 
the external magnetic field. The entanglement shows a quasi-symmetric behavior above and below a 
critical value of the exchange interaction. It becomes more symmetric as the external magnetic field 
increases. The entanglement reaches a large saturation value, close to unity, when the exchange 
interaction is far above or below its critical value and a small one as it closely approaches the 
critical value. Furthermore, we find that the decay rate profile of entanglement is linear when the 
exchange interaction is much higher or lower than the critical value but converts to a power law 
and finally to a Gaussian as the critical value is approached from both directions. The dynamics 
of entanglement is found to be independent of the exchange interaction for isotropically coupled 
quantum dot. 
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I. INTRODUCTION 



Coherence and entanglement lie in the heart of quantum mechanics since its birth and 
are of fundamental interest in modern physics. Particular fields where they play a crucial 
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role are quantum information processing and quantum computi 
is considered as one of the main obstacles toward achieving a practical quantum comput- 
ing system as it acts to randomize the relative phases of the two-state quantum computing 
units (qubits) due to coupling to the environment [5[] . Quantum entanglement is a nonlocal 
correlation between two (or more) quantum systems such that the description of their states 
has to be done with reference to each other even if they are spatially well separated. En- 
tanglement arises naturally when manipulating linear superpositions of quantum states to 
implement the different proposed quantum computing algorithms [6, 7|. Different physical 
systems have been proposed as reliable candidates for the underlying technology of quantum 



There has been 



computing and quantum information processing 
special interest in solid state systems as they can be utilized to build up integrated networks 
that can perform quantum computing algorithms at large scale. Particularly, the semicon- 
ductor quantum dot is considered as one of the most promising candidates for playing the 
role of a qubit 
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171 ]. The main idea is to use the spin S of the valence electron 



on a single quantum dot as a two-state quantum system, which gives rise to a well-defined 
qubit. On the other hand, the strong coupling of the electron spin on the quantum dot to its 
environment stands as a real challenge toward achieving the high coherent control over the 
spin required for information and computational processing. The main mechanisms respon- 
sible for spin decoherence are spin-orbit coupling and spin-nuclear spins (of the surrounding 
lattice) coupling 1?], [3]. Recent experiments show that the spin relaxation time due to 
spin-orbit coupling (tens of ms) is about six orders of magnitude longer than that due to 
coupling to nuclear spins (few ns) 19|, |20|, |21| . This made the electron spin coupling to the 



nuclear spins the target of large number of theoretical 22 , 
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123, 



28 



29 



31 



26J and experimental 



321 ] research works. The dipolar interaction between the nuclei in the 



dot does not conserve the total nuclear spin and as a result the change in the nuclear spin 
configuration happens within T nuc ~ 10 _4 s which is the precession time of a single nuclear 
spin in the local magnetic field of the other spins. This precession time is much longer 
than the decoherence time of the electron spin due to hyperfine coupling in the dot (~ 10~ 6 
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s), therefore we can safely ignore this interaction as well as the fluctuation of the nuclear 
magnetic field when treating the entanglement and decoherence problem of the electron 
spins. In addition to the proposals for using the single quantum dot as a qubit, there has 



been others for using coupled quantum dots as quantum gates [151. |17| . The aim is to find 
a controllable mechanism for forming an entanglement between the two quantum dots (two 
qubits) in such a way to produce a fundamental quantum computing gate such as XOR. In 
addition, we need to be able to coherently manipulate such an entangled state to provide an 
efficient computational process. Such coherent manipulation of entangled states has been 
observed in other systems such as isolated trapped ions [33] and superconducting junctions 



341 ] . Recently an increasing effort to investigate entangled states in coupled quantum dots 



has emerged. The coherent control of a two-electron spin state in a coupled quantum dot 
was achieved experimentally, where the cou plin g mechanism is the Heisenberg exchange in- 
teraction between the electron spins 
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32] . The mixing of the two-electron singlet and 



triplet states due to coupling to the nuclear spins was observed. The induced decoherence 
(leading to the singlet triplet mixing) in the system has been studied theoretically, where 
each one of the two electrons were assumed to localized on its own dot and coupled to a 



different nuclear spin environement 



35 



36 



37]]. There has been proposals for experimental 



scheme to create, detect and control entangled spin states in coupled quantum dots [38 ]. 
Recently, there has been an increasing interest in studying the anisotropy in the exchange 
interaction in coupled quantum dot systems, which is mainly due to the spin-orbit coupling, 
the unavoidable asymmetry of the dots structure and the effect of the external magnetic 



field 
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42J. Furthermore, there has been proposals to utilize this anisotropy, rather 



than removing it, in introducing a set of quantum logic gates to be implemented in quantum 



computing 
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45]. 



In this paper we study the dynamics of entanglement of two electron spins confined 
in a system of two laterally coupled quantum dots with one net electron per dot. The two 
electron spins are coupled through exchange interaction. Modeling the coupling between the 
two electrons on a couple dot system by exchange interaction, including the low barrier case, 
has been studied in detail in Ref. [l^]. We assume that the exchange coupling is anisotropic 
such that the two spins couple in the transverse (z-direction) with a coupling strength that 
is different from the one in the lateral direction (x and y-directions). An external magnetic 
field is applied perpendicularly to the plane of confinement. We assume that the potential 
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barrier between the two dots is low enough that the two electrons are delocalized over the two 
dots and as a result are coupled to a common nuclear environment which splits the two spins 
space into two unmixable subspaces (one contains the singlet state and the other contains the 
three triplets). Therefore, the singlet-triplet mixing can not take place for this system set up 
in contrast to the case studied in Ref . [35) and 36| . Nevertheless, mixing within the triplets 
subspace is possible and in order to study it we assume that the system was initially prepared 
in a maximally entangled triplet spin state with a total antisymmetric wave function, which 



is experimentally an achievable task 3l|, |32j. The two spins are coupled through hyperfine 



interaction to the surrounding unpolarized nuclear spins in the two dots. We study the time 
evolution of the entanglement of the two spins induced by the nuclear spins at different 
strengths of the external magnetic field and exchange interaction, which can be tuned and 
measured experimentally by controlling the different gate voltages on the two dots [32] . We 
found that the entanglement which is initially maximum decays as a result of coupling to 
the nuclear spins within a decoherence time of the order of 10~ 4 . The entanglement decay 
rate and its saturation value vary over wide range and are determined by the value of the 
exchange interaction difference (between the transverse and lateral components) and the 
strength of the external magnetic field. The entanglement exhibits a critical behavior which 
is decided based on the competition between the exchange interaction difference (regardless 
of its sign) and the external magnetic field, it shows a quasi- symmetric behavior above and 
below a critical value of the exchange interaction difference which is determined based on the 
value of the external magnetic field. The behavior becomes more symmetric as the external 
magnetic field increases. Comparable decay rates and saturation values were observed for the 
entanglement on both sides of the critical exchange interaction value. Far from the critical 
exchange value, the entanglement shows high saturation values and slow power law decay 
rates. Close to the critical value, the saturation values become much lower and the decay rate 
shows Gaussian profiles. Our results suggest that for maintaining the entanglement between 
the two spins, the difference between the values of the exchange interaction difference and 
the external magnetic field should be tuned to be large. This paper is organized as follows. 
In sec. II we introduce our model and the calculations of the spin correlation functions. 
In sec. Ill we show the evaluation of entanglement. In sec. IV we present our results, 
discussing and explaining the different important findings we obtained. We close with our 
conclusions in Sec V. 
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II. TWO ELECTRON SPINS COUPLED TO A BATH OF NUCLEAR SPINS 



We consider two coupled quantum dots with a net single valence electron per dot. The 
two electron spins Si and S 2 couple to each other through exchange interaction and couple to 
the same external magnetic field B and to the common unpolarized nuclear spins on the two 
dots {Ii} through hyperfine interaction. The exchange interaction is taken to be anisotropic 
with the coupling strength in the z-direction to be J z and in the x and y directions as Ji. 
The spin-orbit coupling and the dipolar interactions between the nuclei are ignored. The 
system is described by the Hamiltonian 

H = J z S lz S 2z + J\ {Si x S 2x + Si y S 2y ) + gfiB S • B + S • h 2 N7 (1) 

where g is the electron g factor and //# the Bohr magneton. S = Si + S 2 and h 2 N = 
J2i A^I 1 = gA i BH 2 N, where H 2 n is the two dot nuclear magnetic field and the sum is running 
over the entire space. Ai is the coupling constant of S with the ith nucleus and is equal to 
Auol'f/'A^i)! 2 ? where A is a hyperfine constant, t>o(= a 2 a z /N) is the volume of the crystal cell, 
a and a z are the single dot sizes in the lateral and transverse directions, and ipA^i) is the two 
electrons envelope antisymmetric wave function at the nuclear site iv For a typical single 
quantum dot size, the number of nuclei N = 10 6 and their typical nuclear magnetic field 
affecting the electron spin through hyperfine coupling is of the order of ~ A/ (y/Ng fis) 46| . 



which is approximately of magnitude ~ 5 mT [32], with an associated electron precession 
frequency ujn — A/y/N, where lun 3> 1/T nuc , and a ~ 15 nm. The presence of an external 
magnetic field causes Zeeman splitting for the electron spin and in addition has a polarizing 
effect on the nuclear spins. The magnitude of these effects and the associated electron 
precession frequency depends on the strength of the applied field which we consider as a 
varying parameter in this work and is changed over a wide range. Therefore, the precession 
frequency about the external magnetic field can be higher or lower than ujn depending on 
the strength of the external field compared to the nuclear field. We consider some particular 
unpolarized nuclear configuration, represented in terms of the nuclear spin eigenbasis as 
\P Z ), with Jf = ±1/2, where we have considered nuclear spin 1/2 for simplicity. This nuclear 
configurations have the same time evolution as the more general initial tensor product states 
corresponding to arbitrary directions of individual spins [18.]. The Hamiltonian H (Eq. 1) 
is symmetric under exchange of the two electron spins Si and S 2 which reflects the fact 
that in this model the two electrons are considered indistinguishable (delocalized over the 



two dots). This Hamiltonian splits the the two spins space into two subspaces, one of them 
contains the singlet state and the other contains the three triplet states. In this system 
set up (described by the Hamiltonian H), the two subspaces never get mixed. The only 
mechanism by which these two subspaces can get mixed is that each electron spin couples 
to a different nuclear magnetic field (or even a different magnetic field), which is not the 
case here. We assume that the two coupled electron spins are initially (t = 0) in a triplet, 
maximally entangled, state described in the coupled representation by |T ) = | ffJJ- + JJ-ff). 
We are interested in investigating the dynamics of entanglement of this state (dynamics of 
entanglement and its non-monoticity behavior has been studied in many works, for example 
see Ref. [40] and references therein). For that purpose we evaluate the correlators 

Cij(t) = (n|4(*)4'(*)l n )> M = n,y,z ( 2 ) 

where \n) is the total system (two electrons and nuclei) state given by |T )|{/i2}) and Si{t) = 
e ttH SiC~ ltH . The initial state \n) is an eigenstate of the Hamiltonian 

* 2 2 

Ho = J z Siz&2z + -Ji(Si + S2- + S1-S2+) + t z Si Z + h 2 Nz X! ^**> (3) 

^ i=l i=l 

with an eigenenergy e n = — J z /A + Ji/2, where e z = g\isB z ^ consequently we can expand in 
the perturbation (where we use the fact that for a typical nuclear configuration Al/ho n << 1 



which acts as our expansion parameter as can be noticed in the coming calculations 18l. |46||) 



V = "£{S i+ h 2N - + S i „h 2N+ ), (4) 



i=i 



where the total Hamiltonian H = H + V. Using the time evolution operator U(t) 
Texp[— it Jq dt'V(t')], where T is the usual time-ordering operator, we obtain 



Cij(t) = C\f + C iS (t) (5) 



where C\f is the correlator corresponding to V — 0, while CV,(t) is corresponding to the 
corrections due to the perturbation V and are given by 

4 0) = (n\S t S 3 \n) (6) 



and 



C^t) = (nlU^^S^Uit^^S.min) 

fci 

+ E(^l f>t W^Wf>WI^)(A; 2 |[> t (t)S' J (t)[>|n) (7) 
6 



where the intermediate state \h >= |7+i)|{/*}) = | • -Jz 1 = -1/2, • • •}) and 

\k 2 >= = | • -Jz 2 = + 1 / 2 ; ■ ■ •}), while for the singlet state \S Q ) we have 

(Tq\V\S ) = 0. The non- vanishing correlators to the leading order in V read 

S{ 2 = (n\S lz S 2z \n) = -\ + -^{t) + T 2 (t)) , (8) 

S* 2 = (n\S lx S 2x \n) = \, (9) 

S\ 2 = (n\S ly S 2y \n) = \ , (10) 

Ml = {n\S lz \n)= l -(T l (t)-T 2 (t)), (11) 

M 2 z = (n\S 2z \n) = ^(T 1 (t)-T 2 (t)), (12) 

where 

Fi (t) = 4 ]T 1%L S m\u nA t/2) , i = 1, 2 (13) 

where = = A kl /\/2 is the matrix element between initial state |n) = |T )|{- ■ 

= 1/2,- • •}) and the intermediate state \k\) and V nyk2 = A k2 /\/2 is between |n) = 
|T )|{- - -, = -1/2,- • •}> and \k 2 ) and 

VnM = e« - e kl = - (e 2 + (/i 2 ) 2n )] - ^A fcl = cj_ - ^A fcl , (14) 

w„,fe = e n - e fc2 = + (e z + (/i*) 2n )] + = ^+ + T^fe , (15) 

where J = J\ — J z is the exchange interactions difference due to the anisotropy of the coupling 
and (h z ) 2n = (n\h 2 N Z \n) ~ (n\ J2i^ki Il\ n )- We will use h 2n instead of (h z ) 2 n from now on for 
simplicity. As a result of the large value of N, we can replace the sum over ki by integrals 
(E fct F k . = 1/vq JF(t) + 0(1/ N)) and we obtain 

IMO = \[h~ /iCos(cj_ t) - J 2 sin(cj_ *)] , (16) 
r 2 (*) = 4-[ J o - A cos(cj + i) + I 2 sin(cu+ i)] , (17) 



where 

/o = 

^0 



I = — [ dxdydz[A(x,y,z)t] 2 , (18) 
v J 

i = — / dxdydz\A(x,y, z)] 2 cos[A(x,y, z)t/2] , (19) 
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v 

and where 



I 2 = — / dxdydz[A(x,y, z)] 2 sin[A(x,y, z)t/2] , (20) 



A(x,y,z) = Av \?P(x,y,z)\ 2 . (21) 



\ijj(x, y, z)\ 2 is the two electrons envelope antisymmetric wave-function which we construct 
starting from the two single state wave functions (which is the typical wave function used 
to describe the electron state in the confining potential of a quantum dot system 1?], Q) 

&(zi,3/i,zi) = e -(*i-b) 2 /a 2 e -y 2 i/" 2 e -zl/*^ 

2 (* 2 , 1/2, z 2 ) = e -^) 2 l- 2 e -yll- 2 e -*ll< (22) 

which are mixed to obtain an antisymmetric wave function if)(xx, yi, z\\ x 2 , y 2 , z 2 ). Finally 
integrating out the coordinates of one of the two electrons, squaring and normalizing we 
obtain 

|V>(w)| 2 = 3/2 2 ^ h(w/ 2 , e~ 2 y^\- 2 ^e- 2 ^ 2 [cosHAbx/a 2 ) - e^^\ , (23) 

where 26 is the interdot distance. In our calculations we set b = a/2. The integrals Iq, 1% 
and I 2 were evaluated numerically. 



III. EVALUATING THE ENTANGLEMENT OF FORMATION 

In this paper we investigate the time evolution of two-spin entangled state by calculating 
the entanglement of formation E between them. The concept of entanglement of formation 
is related to the amount of entanglement needed to obtain the maximum entanglement state 
p, where p is the density matrix in the two-spin uncoupled representation basis. It was shown 
by WoottersQ that 

E{p)=£{C{p)), (24) 

where the function £ is given by 

S = K l±^*), (25) 
where h(x) = —xlog 2 x — (1 — x)log 2 (l — x) and the concurrence Cis defined as 

C(p) = max{0, Ai - A 2 - A 3 - A 4 }. (26) 



S 



For a general state of two qubits, Aj's are the eigenvalues, in decreasing order, of the Her- 
mitian matrix 

R= ^JVppVp, (27) 

where p is the spin-flipped state of the density matrix p, and defined as 



P= ((Ty® Gy)p*(a y ® Oy). 



(28) 



where p* is the complex conjugate of p. 

Alternatively, the Aj's are the square roots of the eigenvalues of the non-Hermitian pp. 
Since the density matrix p follows from the symmetry properties of the Hamiltonian, the p 
must be real and symmetrical 3] , plus the symmetric property in X and Y direction in our 
system ,we obtain 

( Pi,i ^ 



p 



P2,2 P2,3 
P2,3 P3,3 
^ p 4 ,4 J 



(29) 



The A»'s can be written as 



A Q = h = VPl.lPM, A c = \y/p 2 ,2P3,3 + |P2,3||, A d = \y/p 2 ,2P3,3 ~ |P2,3 



(30) 



Using the definition < A >= Tr(pA), we can express all the matrix elements in the 



density matrix in terms of the different spin-spin correlation functions 



47|, 



50|: 



p ltl = l -M z l + l -M^ + Sl 2 + \ ) (31) 

P2,2 = \m z 1 - 1 -M*-S z 12 + \, (32) 

P3,s = -~M* + - S? 2 + ~ (33) 

P^ = -\Ml- l -M z 2 +Sl 2 + - A , (34) 

P2, 3 = Sf 2 + S\ 2 , (35) 
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Note that if the sign of J is inverted (i.e. J — > — J), uj + — > — oj- and vice versa and as 
a result Ti and T2 exchange their expressions {T\ <-> T2). The only physical quantities that 
are affected by this sign flip are M* and Mf (Eqs. fill I) . (1121) ) which flip sign as well. Despite 
these sign flips the eigenvalues of the matrix R given by Eq. [HU] are not affected because 
under this sign change the density matrix elements p^i <-> and P2.2 P3,3 whereas the 
other elements are invariant. 



IV. RESULTS AND DISCUSSIONS 



In this paper, we focus on the dynamic behavior of entanglement of two coupled elec- 
tron spins under different external magnetic field strengths and at different values of the 
Heisenberg exchange interactions difference J. Interestingly, the entanglement is indepen- 
dent of the sign of the exchange interaction difference J, as was explained in sec. 3, which 
means that regardless of which coupling component is dominant (Ji or J z ) the behavior of 
the entanglement as a function of J stays the same. Firstly, we study the nuclear-induced 
time evolution of entanglement in absence of external magnetic field and with a very small 
anisotropy in the exchange interaction. In fig. [TJ we plot the entanglement of two quantum 
dots versus time with J = A/5000 and e 2 = 0. The initial state of the coupled quantum 
dot is triplet ^=(| ftJ|) + | JJ-ft")), where the two spins are maximally entangled. As shown 
in fig. [1} the entanglement is initially maximum (E=l). Once the interaction between the 
nuclear magnetic field and quantum dots is turned on, the entanglement starts decaying 
with a rapid oscillation and reaches a saturation value of about 0.921 within a decay time 
of the order of 10~ 4 s. As can be noticed the envelope of this oscillating entanglement de- 
creases with a power law profile under this condition before reaching saturation. From the 
process of constructing the density matrix to calculate the entanglement, the entanglement, 
magnetization and spin correlations are directly related. We plot the time evolution of the 
spin correlation function in the z-direction Sf 2 i n the inner panel. The correlation function 
exhibits a very similar behavior to that of the entanglement, it begins with a maximum 
value, -0.25, then decays to a saturation value of about -0.236 within the same decay time. 
Under the first order perturbation theory, the spin correlations in x and y directions are 
unaffected by nuclear spins. This leads to the higher saturation value of entanglement in 
our calculation, which may change at higher orders. Furthermore, because of the initial 
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unpolarized state for the nuclear spins and absence of external magnetic field, the magneti- 
zation of the two quantum dots remains zero. The effect of the external magnetic field with 
a very small anisotropy in exchange coupling between the two spins is shown in fig. [2] where 
e z is set to be A/50 >> h 2n and J = A/5000. The entanglement evolves with a similar 
decaying behavior to what we have observed in fig. [T] except that the decaying amplitude of 
oscillation is much smaller and the saturation value, for both entanglement and correlation, 
is much higher (0.99935 and -0.24989 respectively), which emphasis the important role the 
external magnetic field plays, as expected, significantly enhancing the entanglement between 
the two spins. In fig. [3j we focus on the role played by exchange interaction difference J. 
The dynamics of entanglement is considered at different values of J at two different specific 
external magnetic fields, zero and A/50. In the two panels, J is set to be A/150 and A/130 
in the upper and lower panels respectively. Interestingly, comparing with fig. 1, the satura- 
tion value of the entanglement increases, from 0.92 to 0.997, as we increase J from A/5000 
to A/150 but then decreases back, to 0.93, as we increase J further to A/130. This indicates 
that the larger J does not always help maintaining the entanglement between the two spins, 
i.e. the entanglement is not a monotonic function of J. Obviously, there is an inversion point 
of E between these values of J. The effect of the magnetic field is clear in the three case, 
it suppresses the amplitude of the decaying oscillation and enhances the entanglement as a 
result of its polarizing effect, giving raise to a higher saturation value in the three cases. 
In fact, the entanglement shows very similar behavior at different magnitudes of the exter- 
nal magnetic field and exchange interaction difference, however the asymptotic saturation 
values and decay rates vary significantly based on the values of J and e z . In principle, each 
electron spin is in a precession motion about two fields, magnetic which is uniform and 
the non-uniform nuclear field. This precession motion is the reason behind the oscillatory 
behavior of the entanglement, whereas coupling to the non- uniform nuclear field is respon- 
sible for the decay where the spin energy is dissipated to the nuclear environment (the same 
behavior was discussed in detail in Ref. [18] for a single spin on a single dot). Therefore, 
it is of great interest to study the mutual effect of the exchange interaction difference and 
the external magnetic field on the entanglement. The phase diagram of the system is shown 
in fig. H] where The saturation value of entanglement is plotted versus J at different values 
of the external magnetic field. As the entanglement mostly reaches its saturation value at 
t ~ 10 _4 s, we use the value of E at t — 10 -4 to approximately represent the saturation 



11 



value at each J. There are three pairs of curves with each pair corresponding to a specific 
value of the external magnetic field given by, from left to right, A/50, A/500, 0. The pair 
of curves corresponding to e z = A/ 50 are replotted in the inner panel for clarity. It is very 
interesting to note that for each e z the behavior of the entanglement is almost symmetric 
about a critical value of J, we call it J c . On each pair, the entanglement saturation value 
is initially unity or very close to unity , depending on the external magnetic field strength, 
then starts to decrease as J increases. It decays exponentially as J approaches J c but once J 
exceeds J c it raises up exponentially again reaching a value of unity or close. The decay rate 
of the entanglement saturation value for J > J c is more rapid than J < J c but they become 
more symmetric as e z increases. Surprisingly, this phase diagram shows a very peculiar role 
played by the exchange interaction difference, namely that small values (J << J c ) and great 
values as well (J >> J c ) enhances the entanglement between the two spins, while as J — > J c 
it is reduced significantly. 

Unfortunately, we can not investigate the entanglement behavior when J is very close or 
equal to J c because the system diverges under this condition. The reason for this divergence 
is the degeneracy (electronic levels crossing) that takes place between the initial state and 
the intermediate states which causes a break down of the first-order perturbation theory. 
This degeneracy is manifested in fig. [5J where the energy difference ui n! ki between the initial 
state |n) and the first intermediate state \ki) and the difference cu n ,fc2 between \n) and \k 2 ) 
(given by Eqs. [T4l and [TBI) are plotted versus the exchange interaction difference J for chosen 
h>2n = A/ 500 and e z = A/ 50. As can be seen in fig. [5j uj n ^\ vanishes at a particular value 
of J which turns out to be the same as J c observed in the inner panel in fig. H] where 
the same values of parameters are used in that case. Since io n ^\ appears squared in the 
denominator of the first order perturbation term (Eq. [T3l) . it is responsible for the divergence 
when it vanishes. Also, the symmetric behavior of the entanglement about J c is due to the 
linear dependence of the magnitude of iO n ,ki on J as h increases or decreases away from J c . 
Similarly, co> n ,fc2 may cause the same effect if negative values of J were considered. The value 
of the critical exchange interaction J c corresponding to each e z can be easily deduced from 
the phase diagram to be equal to 2(h 2n + e z ) coinciding with the observations of fig. [51 which 
emphasis that the external magnetic field does not only enhance the entanglement but also 
dictates the value of J c . 

After the important observations from the phase diagram (fig. Hj), in particular concerning 
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the entanglement critical behavior upon varying J, it is interesting to further investigate the 
effect of J on the entanglement above, below and close to J c at a fixed value of the external 
magnetic field. Therefore, the dynamics of entanglement at different exchange interaction 
values below J c in fig. Eland above it in fig. [7] is examined. Where e z is set to be A/100 
and consequently J c turns out to be A/41. 7. In fig. Owe consider the exchange interaction 
values A/45, A/50 and A/55 (< J c ) whereas in fig. [7J A/40, A/37 and A/35 (> J c ). As we 
can see in the two graphs, when the value of J is far from J c , either higher or lower, the 
entanglement oscillation is small and the decay rate is slow reaching high saturation values 
that approach unity as J becomes much smaller or much higher than J c . As J gets very close 
to J c the entanglement exhibits very rapid oscillation with very large initial amplitude and 
decays very rapidly as well reaching a low saturation value, 0.74, for J < J c and very low one, 
0.22, for J > J c . The contrast between those values gets smaller and smaller with increasing 
the magnetic field as can be noticed from the phase diagram (fig. HJ). To compare the decay 
rates of entanglement at different J values, the corresponding entanglement decay profiles 
are plotted in the inner panels in figs. [6] and UJ As shown, the decay profiles are linear when 
J is much higher (or lower) than J c converting to a power-law and then to a Gaussian as J 
gets closer to J c from above or below. The critical behavior of the entanglement based on 
the strength of J points out that there is a competition between the effects of the external 
magnetic field, which enhances the entanglement through polarizing the two spins in its 
direction, and the exchange interaction difference, which enhances the entanglement on the 
other hand by favoring anti-ferromagnetic correlation. Therefore when one of the two effects 
is dominant, i.e. J « e z (or >> e z ), there will be a large net correlation one way or the 
other leading to strong entanglement. Tuning the two effects to be comparable reduces the 
net correlation significantly and forces the entanglement to attenuate to a low value. 



V. CONCLUSIONS 



In Summary, we study the time evolution of entanglement of two electron spins in 
anisotropically coupled quantum dot under coupling with nuclear spins in an applied ex- 
ternal magnetic field in the triplet (S = 1) subspace. Our main observations is that the dy- 
namics of entanglement exhibits a critical behavior determined by the competition between 
the value of the exchange interaction difference between the two spins and the strength of the 
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external applied magnetic field. When one of them dominates, the entanglement saturates 
to a very large value close to the maximum. On the other hand, when they are comparable 
the entanglement attenuates and may reach very low saturation values. We found that for 
each fixed value of the external magnetic field there is a corresponding critical value of the 
exchange interaction difference where the entanglement between the two spins as a function 
of the exchange interaction exhibits symmetrical behavior about it. When we studied the 
time evolution of the entanglement at values of the exchange interaction difference above and 
below the critical value, we observed that the entanglement decay profiles are linear when 
the exchange interaction value is much higher or lower than the critical value and saturates 
to large values but then converts to a power law and finally to a Gaussian as it approaches 
the critical value from above or below reaching low saturation values. The system diverges 
when the exchange interaction value becomes very close or equal to the critical value which 
might be a sign of a critical behavior taking place in the environment consisting of the 
nuclear spins. For an isotropic coupling the entanglement becomes absolutely independent 
of the Heisenberg exchange interaction. Our observations suggest that tunning the applied 
external magnetic field and (or) the exchange interaction difference between the two spins 
in such a way that they differ significantly in value may sustain the correlation between the 
two spins in the coupled quantum dot and lead to a strong entanglement between them. 
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x 10 



FIG. 1: Time evolution of entanglement and spin correlation in z-direction in absence of external 
magnetic field and for very small anisotropy J = ^4/5000. The nuclear magnetic field hm = ^4/500, 
where A is the hyperfine constant. 
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FIG. 2: Time evolution of entanglement and spin correlation in the z-direction for very small 
anisotropy J = ^4/5000 with a large external magnetic field e z = A/50 » hm- 
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FIG. 3: Dynamics of entanglement at different exchange interaction differences ^4/150, and ^4/130 
with the external magnetic field e 2 set at and A/ 50. 
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FIG. 4: The phase diagram of entanglement. The saturation value of the entanglement is plotted vs. 
the exchange interaction difference J at different values of the external magnetic field e z = 0, A/500, 
and A/50. The case corresponding to e z = A/50 is replotted in the inner panel to demonstrate the 
symmetric behavior of E at high external magnetic field. 
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FIG. 5: The energy differences w n> fci (solid line) and u n ^2 (dashed line) are plotted versus the 
exchange interaction difference J for chosen h,2 n = ^4/500 and t z = A/ 50. 
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FIG. 6: Dynamics of entanglement with the external magnetic field fixed at e z = ^4/100 while the 
exchange interaction difference J is set at values below the critical value (J c = ^4/41.7), J = A/45, 
A/ 50 and A/ 55. The inner panel shows the corresponding entanglement decay profiles. 
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FIG. 7: Dynamics of entanglement with the external magnetic field fixed at e z = ^4/100 while the 
exchange interaction difference J is set at values above the critical value (J c = A/41.7), J = A/35, 
A/ 37 and ^4/40. The inner panel shows the corresponding entanglement decay profiles. 
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